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Quantum tunneling of the superconducting order parameter gives rise to the phase slippage 
process which controls the resistance of ultra-thin superconducting wires at sufficiently low temper- 
atures. If the quantum phase slip rate is high, superconductivity is completely destroyed by quan- 
tum fluctuations and the wire resistance never decreases below its normal state value. We present 
a detailed microscopic theory of quantum phase slips in homogeneous superconducting nanowires. 
Focusing our attention on relatively short wires we evaluate the quantum tunneling rate for phase 
slips, both the quasiclassical exponent and the pre-exponential factor. In very thin and dirty metallic 
wires the effect is shown to be clearly observable even at T — > 0. Our results are fully consistent with 
recent experimental findings [A. Bezryadin, C.N. Lau, and M. Tinkham, Nature 404, 971 (2000)] 
which provide direct evidence for the effect of quantum phase slips. 



I. INTRODUCTION 



It is well established that superconducting fluctuations play a very important role in reduced dimension. Above 
the critical temperature Tc such fluctuations yield an enhanced -conductivityEJ. Below Tc fluctuations are known to 
destroy the long-range order in low dimensional superconductorsEI. Does the latter result mean that the resistance of 
such superconductors always remains finite (or even infinite), or can it drop to zero under certain conditions? 

It was first pointed out by Littled that quasi-one-dimcnsional wires made of a superconducting material can acquire 
a finite resistance below Tc of a bulk material due to the mechanism of thermally activated phase slips (TAPS). This 
TAPS process corresponds to local destruction of superconductivity by thermal fluctuations. Superconducting phase 
tp(t) can flip by 2tt across those points of the wire where the order parameter is (temporarily) destroyed. According 
to the Josephson relation V — ip/2e (here and below we set h = 1) such phase slips cause a nonzero voltage drop 
and, hence, dissipative currents inside the wire. A theory of this TAPS phenomenon was developed in Refs. 0M This 
theory yields a natural result, that the TAPS probability and, hence, resistance of a superconducting wire R below 
Tc are determined by the activation exponent 

R(T) « exp(-£//T), U ~ N ° A ^ T) s g(r), (1) 

where U (T) is the effective potential barrier for TAPS determined simply as the superconducting condensation energy 
(No is the metallic density of states at the Fermi energy and A(T) is the BCS order parameter) for a part of the wire 
of a volume s£ where superconductivity is destroyed by thermal fluctuations (s is the wire cross section and £(T) is 
the superconducting coherence length). At temperatures very close to Tc ew. (|l|) yields appreciable resistivity which 
was indeed detected experimentallyQ. Close to Tc the experimental resultsa fully confirm the activation behavior of 
R(T) predicted in eq. (|l|). However, as the temperature is lowered further n lielow Tc the number of TAPS decreases 
exponentially and no measurable wire resistance is predicted by the theoryem except in the immediate vicinity of the 
critical temperature. 

Experimentstj were done on small diameter whiskers and thin film samples of typical diameters ~ 5000 angstrom. 
Recent progress in nanolithographic technique allowed to fabricate samples with much smaller diameters down to 
~ 10 nm. In such systems one can consider a possibility for phase slips to be created not oidy due to thermal but also 
due to quantum fluctuations of a superconducting order parameter. Mooij and coworkers!] discussed this possibility 
and attempted tOrObserve quantum phase slips (QPS) experimentally. 

Later GiordanoH performed experiments which clearly demonstrated a notable resistivity of ultra-thin supercon- 
ducting wires far below Tc- There observations could not be adequately interpreted within the TAPS theory and were 
attributed to QPS. Other groups also reported noticeable deviations from the TAPS prediction in thin (quasi-) ID 
wires tM " 

First theoretical studies of the QPS effectslilliiia were performed within a simple approach based on the time- 
dependent Ginzburg-Landau (TDGL) equations. Later in Refs. [lJlPi a microscopic theory of QPS processes was 
developed with the aid of the imaginary time effective action techniqueH which properly accounts for non-equilibrium , 
dissipative and electromagnetic effects during a QPS event. One of the main .conclusions reached in Refs. |l4|,[l5]is that 
the QPS probability is considerably larger than it was predicted previouslyE3. For ultra-thin superconducting wires 
with sufficiently many impurities and with diameters in the 10 nm range this probability can already be large enough 



1 



to yield experimentally observable phenomena. Also, further interesting effects including quantum phase transitions 
caused by interactions between quantum phase slips were discussedoEj. r, 

In spite of all these developments an unambiguous interpretation of the resultsB in terms of QPS could still be 
questioned because of possible granularity of the samples used in these experiments. If that was indeed the case, QPS 
could easily be created insidej-weak links connecting neighboring grains. A.lso in this case superconducting fluctuations 
play a very important roleEDtJ, however - in contrast to the QPS scenarioEStJ - the superconducting order parameter 
needs not to be destroyed during the QPS .event. 

Recently, Bezryadin, Lau and Tinkhamo developed a new technology which allowed them. to fabricate essentially 
uniform superconducting wires with thicknesses down to 3-^5 nm. According to our theoryi^lilthe QPS effects should 
be sufficiently large in such systems to be observed in experiments. And indeed, the authorsEj observed that several 
wires showed no sign of superconductivity even at temperatures well below the bulk critical temperature. Moreover, 
at lower temperatures their resistance was found to increase with decreasing temperature, i.e. these samples could 
even turn insulating at T — > 0. The authoK^a also argued that thep, experimental data can be interpreted in terms 
of a quantum dissipative phase transition 2 -*^ which was predictedta also for ultra-thin superconducting wires in a 
certain parameter range. r-. 

The resultsE3 are qualitatively consistent with previous experimental findingsO. Both experimental works support 
our general understanding of the role of QPS processes in mesoscopic superconducting wires and call for more detailed 
theoretical studies of the QPS effects. In Refs. an importance of collective modesE3 and QPS interaction effects 

was mainly emphasized. These are particularly important for long wires. On the other hand, for relatively short 
wires interaction between different phase slips - at least its spatially dependent part - should not play any significant 
role. Let us now recall that the wires studied in the experimentsEHl are not only considerably.-thmner but also much 
shorter than those investigated by GiordancQ. To give some numbers, the length of the wiresO was typically 40 -j- 50 
/zm whereaSpthe wiresca were only 0.1 -j- 0.2 /xm long. At the same time, the superconducting coherence length in the 
experiments 2 !] was even shorter, £ ~ 7-j- 8-prn, i.e. such samples can still be considered as quasi-lD superconductors. 

Motivated by the experimental findingscj, in this paper we will present a detailed microscopic investigation of single 
quantum phase slips. We will focus our attention Wrafi accurate evaluation of the QPS tunneling rate rather than on 
the interaction effects between different phase slipsllirll3. We will go beyond the exponential accuracy and also evaluate 
a pre-exponential function in the expression for the QPS rate. We will then use our results for a direct quantitative 
comparison with the experimental results 2 ^. 

The structure of the paper is as follows. In Section 2 we will formulate a simple derivation of the effective action 
for our problem with an emphasis put on the Ward identities. In Section 3 we will make use of our general results 
and derive the action for a special case of ultra-thin superconducting wires. We also evaluate the QPS rate within 
the exponential accuracy. Section 4 is devoted to an estimate of the pre-exponent for this rate. Comparison with 
experiments and brief conclusions are presented in Section 5. Some further technical details are diverted to Appendix. 



II. THE MODEL AND EFFECTIVE ACTION 



The starting point for our analysis is a model Hamiltonian that includes a short range attractive BCS and a long 



range repulsive Coulomb interaction. The idea is to integrate out the electronic degrees-' 
partition function, so that we are left with an effective theory in terms of collective fields 



sdom on the level of the 
The partition function 



Z is conveniently expressed as a path integral over the anti-commuting electronic fields if), tp and the commuting gauge 
fields V and A, with Euclidean action 



S = I dx(^p a [d r 

-^ipfipliplipt +ienV 



ieV + Z{V- -A)h/>„ 
c 



[E 2 +B 2 ]/8tt 



(2) 



Here £(V) = — V 2 /2m — u + U(x) describes a single conduction band with quadratic dispersion and also includes an 
arbitrary impurity potential, A is the BCS coupling constant, a =|, J, is the spin index, and en denotes the background 
charge density of the ions. In our notation dx denotes c£ 3 xg?t and we use units in which h and fcs are set equal to 
unity. The field strengths are functions of the gauge fields through E = —W + (l/c)d T A and B = V x A in the 
usual way for the imaginary time formulation. 

We use a Hubbard-Stratonovich transformation to decouple the BCS interaction term and to introduce the super- 
conducting order parameter field A = Ae' v 

V 2 Ae-if dxA2 



exp A / dx^ipiipi^ 
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x /*X' 2 Ae-/ d;E (* A2+A *^i+ A >^T) ; ( 3 ) 



where the first factor is for normalization and will not be important in the following. As a result, the partition function 
now reads 

Z = Jv 2 Aj V 3 A J VVVHei- 30 -!^ 5 "*) , (4) 

//p2 I d2 a 2 \ 

dxl ±- + i en V+ — \ , 

where the Nambu spinor notation for the electronic fields and the matrix Green function in Nambu space 



C-i= [ ^r-iey + e(V-|A) A \ 

y A* 9 T + ie y -£(V + fA) ) " U 

has been introduced. After the Gaussian integral over the electronic degrees of freedom, we are left with the final 
effective action 

S eS = -Tr]n(}- 1 +S [V,A,A]. (6) 

Here the trace Tr denotes both a matrix trace in Nambu space and a trace over internal coordinates or momenta and 
frequencies. In the following "tr" is used to denote a trace over internal coordinates only. 

The gauge invariance of the theory enables us to rewrite the action ([|) in a different form, which is more convenient 
for us, 

5^ = -Trlng- 1 +5 [V,A,A], (7) 

where 



C-i = ( 9 T +e(V)-ie$+^-|{V,v s } A 

V A S T -e(V)+ie$-^-i{V,v s } I ' 



and we have introduced the gauge invariant linear combinations of the electromagnetic potentials and the phase of 
the order parameter 

* = V-L v . = -L(v<p-^Ay (9) 



2e 2m \ c 

The curly brackets {A, B} denote an anti-commutator 



A. Perturbation theory 

The action (^) cannot be evaluated exactly. Here we will perform a perturbative expansion in $ and v s . We will 
keep the terms up to the second order in these values. This perturbation theory is sufficient for nearly all practical 
purposes, because nonlinear electromagnetic effects (described by higher order terms) are known to be usually very 
small in the systems in question. Our general derivation holds for an arbitrary concentration and distribution of 
impurities as well as for arbitrary fluctuations of the order parameter field in space and time. 

We split the inverse Green function ([|) into two parts 

g -x_(d T +m a \ 



and 



1 I o ,< -I- " • ' 
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The logarithm in the equation (Q) can now be expanded in powers of Q l 1 and we get 

Trln^- 1 = Trln^ 1 + Tr(So£q 1 ) - ^Tv(g g^) 2 . (12) 
The Green function Go has the form 

'%%). (13) 

In eq. ([l^) we used the fact that the non-diagonal component A in the matrix Q^ 1 is real. As a result we have F = F, 
F(xi,X2) — F(x2,xi) and G{x\,X2) = —G(x2,xi). 

B. Ward identities 

The Green function Go satisfies an important identity, which is easy to check: 

^X-X^^g-^g}^, (14) 

where x ls an arbitrary function of time and space, and 03 is one of the Pauli matrices. Multiplying this matrix identity 
by Q from the left and from the right side and taking the diagonal components of the resulting matrix equation we 
get two identities: 

xG-g x = g(x-{^})g + f( x + { v^})f, 

X G-Gx = f(x- {V, F + G (x + {V, G. (15) 

Below we will use these identities in order to decouple the effective action of the BCS superconductor and to reduce 
it to a transparent and convenient form. It is important to emphasize again that these identities are valid for any 
impurity distribution and for any time and spatial dependence of the order parameter field. It is also worth mentioning 
that the Ward identity ( pT} ) is not the result of the gauge invariance of the theory. It remains valid even for uncharged 
particles. 

The Ward identity related to the gauge invariance of our theory has a different form: 

Go X <r 3 X - X^Go 1 = ?r°3 - {V, ^} - 2ta 2 A x . (16) 
or 2m 

We will use this identity to transform the first order correction to the action. It is interesting, that in the absence of 
superconductivity the identities ( |l4| ) and (|l6|) are equivalent because the inverse Green function commutes with 03 in 
this case. For superconductors, however, these two identities are different. 



C. First order 



The first order correction to the effective action is 



Si = -TtiQoGi 1 ) = -tr 



771V" 



j e $ ( G _G)--{V,v s }(G + G) 



(17) 



With the aid of the Ward identity (jlfj) it easy to show that the phase of the order parameter drops out from the first 
order terms in the electromagnetic fields. The action Si can therefore be rewritten as 



Si = -tr(mv^G) 



dx 



1 



A]T/+-j e [A]A . 



(18) 



We note that in general the electron density 7i e [A] and the current density j e [A] explicitly depend on the absolute 
value of the order parameter. 
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D. Second order 



It is convenient to introduce the following notations: 

9 = 2e<I>, £ = {V,v s }. 
In terms of these new variables the second order correction to the action reads 



(19) 



1. 



<S 2 = 2 Tr (£o£i 



-X\2 



1 

-~tr 



G9G9 + G9G9 - 2F9F9- 
GCGC + GCGC + 2FCFC + 2G9GC - 2G9GC 



(20) 



Here we used the properties of the Green function (|13|). The form (|20|) of the second order correction is not quite 
convenient, because it contains the 0C terms. In order to separate 9 and C we use the Ward identities (|l5|). We write 

G9G = 9G-G9 + G{V, — }G -F\9 + {V. — } ) F, 

2m \ 2m J 



GOG = 8G-G9 + G{V, - F (o - { V, ^}) F. 

Inserting these expressions into (pp| ) after some simple transformations we rewrite the second order contribution as 
follows 



S 2 = -tr(G(v.V0)) - I 



GKGK + GKGK - GMGM + GMGM - G9G9 



Here we have introduced 



-G9G9 - 2FKFK - 2F9F9 + 2FMFM + AFCFC 



^ = {V,u}, M = {V,^}, 



(21) 



The values 9 and v s are now almost decoupled. The terms containing both these values were transformed into the 
terms containing the linear combination of these values u which does not depend on the phase of the order parameter 
field. The action @ can be simplified further. We rewrite the first term of (|TJ) as follows 



-tr(G(v s V0)) = tr 



Am 



in 2 I G 



Again we decouple 9 and v s . Making use of the identities ( |l5| ) yet a couple of times we arrive at the final expression 
for the second order contribution to the effective action: 



S 2 = tr{mx 2 s G) - tr(mu 2 G) - -tr(G{V, u}G{V, u}) 
itr(F{V, u}F{V, u}) + l -tr{F9F9) ~ itr(F{V, v s }F{V, v s }). 



(22) 
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E. Resulting action 

Combining all contributions, we get the final resultEl 

S = 5 s [A,$,v s ] + SW[A, V, A] + S em [E, B] , 

where 



(23) 



S s = J dx (^pj - THnSo X [A] + Trln^A = 0] 



^tv(F6F9) - itr(F{V,v s }F{V,v s }), 



(24) 



S N = J dx \-ie(n e [A] - n)V - ij e [A]A- 



mu 



5 [A] 



~tr(G{V, u}G{V, u}) + Jtr(F{V, u}F{V, u}), 



dx 



E 2 + B 2 

8^ 



(25) 
(26) 



III. EFFECTIVE ACTION FOR ULTRA-THIN WIRES 



A. Averaging over the electromagnetic field 



The above expressions are complicated and in general can hardly be evaluated in a closed form. In this section 
we will focus our attention specifically on the case of quasi-one-dimensional superconducting wires and calculate the 
effective action performing several approximations. We will argue that our procedure allows to evaluate the QPS 
action up to a numerical prefactor of order one. 

If one assumes that deviations of the amplitude of the order parameter field from its equilibrium value are relatively 
small, the above effective action can be significantly simplified. We expand the general effective action ( p3| ) in powers 
of SA(x,t) — A(x,t) — A (hece A = Abcs) up to the second order terms. The next step is to average over the 
random potential of impuritiesc°l. After that the effective action becomes translationally invariant both in space and 
in time. Performing the Fourier transformation we obtain 




(27) 



Here L and C are respectively the inductance times unit length and the capacitance per unit length of the wire. The 
functions \e, XJi Xl and \A, which depend both on the frequencies and the wave vectors, are expressed in terms of 
the averaged products of the Green functions appearing in the eqs. (^3), ( p5| ) (see Ref. |l6| for more details). These 
functions can be evaluated analytically for most limiting cases. For the sake of completeness some explicit expressions 
are presented in Appendix. 

The voltage V and the vector potential A enter the action in a quadratic form and, hence, can be integrated out 
exactly. After that the effective action will only depend on ip and 5 A. We find 



b ~ 2 



dtodq 



{.F^I^ + x^AI 2 }, 



(28) 



where 



G 



, ,2 j_ Xl „1 



XJ, 



X 



XE 



Cto 2 



XjXl 
4m 2 



Cu 2 



13 + xj 



XEQ 2 ) (±+Xeu 2 
The electromagnetic potentials are expressed as follows: 



tXl 



(29) 



V 



A 



XJ[7T+ Xeu 



-XL 



tXeXlQ 



XJ + XEQ Z ) 

rn?XL (~ 



' 1 

. sL 



Xe^ 



xj + xeq 



n*XL) 



-IbJ 



T+XJ + XEq 2 ) + XE^ 2 + ^ X l) 



2e 



2e 



(30) 
(31) 



In most of the situations the wire inductance is not important and can be neglected. Therefore here and below we 
put L = 0. Then we get 



todq j(^ 2 + ^q 2 )^ +XE<?)+_%&£ 2 



(2n) 2 



XJ + X£<? 



Xa\SA\ 



and 



V = 
A = 



XJ 



xj + XEq z 



2e 



0. 



(32) 

(33) 
(34) 



Let us note that the Josephson relation V — <p/2e is in general not satisfied. According to eq. ( p3|) this relation may 
approximately hold only in the limit XJ 3> C/s + XEq 2 ■ Making use of the results presented in Appendix one easily 
observes that in a practically important limit of small elastic mean free paths / the latter condition is obeyed only at 
low frequencies and wave vectors u/A <C 1 and Dq 2 / A <§C 1, where D = vpl/3 is the diffusion constant. 

Let us now perform yet one more approximation and expand the action in powers of uj and q 2 . Keeping the terms 
of the order q A and cu 2 q 2 we find 



S 



dudq 
+2N 



C 



iraAq 



-aDq* 



7T(T 

8A 



2„2 



oj q 



) 
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2e 



12A 2 



irDq 2 
8A 



l*A| 



(35) 



The term oc turns out to be equal to zero. In Pq ) we introduced the normal state conductance of the wire 
a = 2e 2 N D. At even smaller wave vectors, Dq 2 /2A <C 2C/-Ke 2 N s <C 1, we get 



b ~2 



dtodq 



(Cuj 2 + iraAsq 2 



2c 



sXa\SA\ 



(36) 



Here we have assumed C/2e 2 Nos -C 1. This inequality is usually well satisfied for sufficiently good metals, perhaps 
except for the case of some specially chosen substrates. The form of the action suggests the existence of the plasma 
modes which can propagate along the wire. These are the so-called Mooij-Schon modesES, the velocity of which is 
given by the equation: 



Co 



iraAs 



C 



(37) 



B. QPS action 

One can showQ that for very long wires the action (^6|) yields a QPS solution described by a simple formula 
ip(x, t) = — arctan(z/coT). The long time behavior of this solution results in the logarithmic interaction between two 
phase slips (xi,t%) and (x2, t^): 
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Sin 



■In 



(zi -x 2 ) 2 +cg(r 1 -t 2 ) ; 
£ 2 



(38) 



where 



' sC 



(39) 



Here a ~ 1/137 is the fine structure constant. In short wires, however, the above logarithmic interaction (|3§|) does 
not play an important role and can be essentially neglected. 

Let us estimate the contribution of a single phase slip to the effective action. First we rewrite the action ( |35| ) in 
the space-time domain dropping the unimportant term oc (Dq 2 ) 2 : 



s = U dxdT \Ss 



dip 



■no A 



dp 
dx 



TT(J 



32e 2 A \dxdr 



sN dxdrl 5 A 2 



i fas Ay ttD fdSAY 

YlAtXjh) + 8A\dx~) 



(40) 



Then we assume that the absolute value of the order parameter is equal to zero at a time r = and at a point 
x = 0. The size of the QPS core is denoted as xq, and its time duration is To. The amplitude of the fluctuating part 
of the order parameter field \8A(x, r)| can be approximately expressed as follows: 



\5A(x,t)\ = Acxp{-x 2 /2xl -t 2 /2t 2 ). 



(41) 



The QPS phase dependence on x and r should satisfy several requirements. In a short wire and outside the QPS 
core the phase tp should not depend on the spatial coordinate in the zero current bias limit. On top of that, at x = 
and r = the phase should flip in a way to provide the change of the net phase difference across the wire by 2ir. For 
concreteness, let us present two different trial functions which obey the above requirements. For instance, one may 
choose 



<p(x,r) 



2 cosh(r/ro) 



tanh 



xo tanh(r/ro) / ' 



(42) 



ip(x,r) 



— tanh ( 

2 \x t 



(43) 



Similar other trial functions can also be considered. 

Substituting the trial functions ([i"l|), ( |4^ ) (or (41), J43])) into the action j4C| ) one arrives at the expression 



S(xq,t ) 



Ol- 



C 



a 2 sN Q 



a 3 sN DA 



xo 



sN oD 1 A7 a2 C 

■ a4 — t h a 5 sN Q A x t + a 6 

A x T 



9 ! 

e 2 T 



(44) 



where a,j are numerical factors of order one which depend on the precise form of the trial functions, C = CX is the 
total capacitance of the wire and X is the wire length. Note that fictitious divergences emerging from a singular 
behavior of the functions (^J) , d43| ) at x — xq and t — tq are eliminated since the order parameter vanishes inside the 
QPS core. 

Let us first disregard capacitive effects neglecting the last term in eq. (|44|). Minimizing the remaining action with 
respect to the core parameters xq and To and making use of the inequality C/e 2 Nos <C 1, we obtain 



•I'D 



/ 0304 



1/4 



A' To 



a 2 a4 \ 1 



\ a 2a5 J V A V a 3 a 5 

These values provide the minimum for the QPS action, and we find 



(45) 
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Sqps = 2(V^3 + Va^)N sv r DA. (46) 
One can also express eq. ( ff6| ) in the form convenient for further comparison with experiments: 



S QPS = A^. (47) 



Here A = 2(^/0203 + ^/ 0,4125) /tt, R is the total wire resistance, R q = 7r?i/2e 2 = 6.453 kf2 is the resistance quantum 
and £ = \jDj A is the superconducting coherence length. 

As it was already pointed out, the results ( |45| ) and ([l6]) hold provided the capacitive effects are small. This is the 
case for relatively short wires 



In the opposite limit the same minimization procedure of the action ( f4~i| ) yields 



X«t^. (48) 



J'O 



£, Ar ~ y/XC/^Nos » 1. (49) 

The QPS action again takes the form @ with yl - XC / £,e 2 N s. 

For the sake of clarity, let us summarize the approximations performed in this section. As a first step, we expanded 
the action derived in Section 2 up to the second order in 8A[x,t) = A(x,t) — A. Obviously this approximation is 
sufficient everywhere except inside the QPS core where A(x, r) is small. In these space- and time-restricted regions 
one can expand already in A(x,t) again arriving at eq. ( |27|) with 5A(x,t) —> A(x,t) and with all the x-functions 
defined in Appendix with A = 0. Both expansions match smoothly at the scale of the core size xq ~ £, To ~ 1/A. 
Hence, the approximation ( |27| ) is sufficient to obtain the correct QPS action, perhaps up to a numerical prefactor of 
order one. 

In order to simplify our analysis further, in (35) we expanded ( p7j ) in powers of uj/ A and Dq 2 /A. Again this 



approximation is sufficient within the same accuracy. Indeed, one can - even without performing this expansion 
- substitute the trial functions (|4l|), ( f42| ) (or ([fl|), (f43|)) directly into the action (|32|). If the capacitive effects are 
neglected (ff8|), the resulting QPS action can be represented as a function of the dimensionless parameters £o/£ and 
Atq only. Making use of the general expressions for the x-functions collected in Appendix and minimizing the QPS 
action with respect to Xq and r one again arrives at the result (^) with A ~ 1. If the inequality dig) is violated, the 
accuracy of our expansion in powers of u)/A may only become better (cf. eq. ((49|)). 

Finally, the particular choice of the trial functions (e.g. ([|l]) and ([l2|)) describing the QPS event also appears not to 
play any significant role as long as these trial functions obey the general requirements formulated above. In addition 
to ©, d|) and (|f) we have used several other trial QPS functions. In all cases we have obtained A within the 
interval A 1=3 0.8 -j- 2.5. In a way our method can be regarded as a variational procedure. Therefore, even though the 



exact value of a numerical prefactor A in eq. (47) cannot be established within our approach, we do not expect A to 
deviate substantially from the above values. 

Our last remark in this section concerns the role of dissipation. From the form of the result (^6|) one could naively 
assume that the correct QPS action could be guessed, e.g., from a simple TDGL-based approach (or, alternatively, only 
from the "condensation energy" term proportional to xa) without taking into account dissipative and electromagnetic 
effects. Indeed, minimization of the contribution ~ |<5A| 2 (the last three terms in eq. (|40|)) is formally sufficient to 
arrive at the correct estimate Sqps ~ N^sy/ DA. It is obvious, on the other hand, that not only the amplitude but also 
the phase fluctuations of the order parameter are important during the QPS event. If the latter fluctuations are taken 
into account without including dissipative effects (this would correspond to formally setting a — > in eq. ([40|)) the 
estimate for the QPS action Sqps ~ fJ- would follow immediately. This result would be parametrically different from 
eq. ([46|). Note, however, that within our model the dissipative effects can be ignored only for C/e 2 Nos > 1. Usually 
the latter condition cannot be satisfied for metallic systems, perhaps except for some specially chosen substrates. In 
the opposite - more realistic - limit C/e 2 Nos 1 dissipation plays a dominant role during the phase slip event, 
and the correct QPS action cannot be obtained without an adequate microscopic description of dissipative currents 
flowing inside the wire. 

IV. PRE-EXPONENT 

The above results allow to estimate the exponential suppression of QPS in ultra-thin superconducting wires de- 
pending on thickness, impurity concentration and other parameters. These results, however, are not yet sufficient to 
evaluate the whole QPS rate which has the form 
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Iqps = Bexp(-S QPS ). 



(50) 



The task at hand is to provide a reliable estimate for the pre-exponential factor B in eq. ([50j). A general strategy to 
be used for this purpose is well knowncB. One can start, e.g., from the expression for the grand partition function of 
the wire 

Z = J VAVipexp(-S) (51) 

and evaluate this path integral within the saddle point approximation. The least action paths 

SS/\5A\ = 0, 6S/6<p = (52) 

determine all possible QPS configurations. Integrating over small fluctuations around all QPS trajectories one repre- 
sents the grand partition function in terms of infinite series (each term in such series corresponds to one particular 
QPS saddle point). Then - at least if interaction between different quantum phase slips is small and can be neglected 
- one can easily sum these series and represent the final result in the form of the exponent 

Z = exp(-F/T), (53) 

where a formal expression for the free energy F reads 

^ JVSYexpj-S^iSY}) lQPS ( 

T JV6YeM^S [SY}) ^~ S Q^) = F ° ~ ~ ^ 

Here F is the free energy without quantum phase slips, SY — (5 A, Sip) denote the fluctuations of relevant coordinates 
(fields), ^ 2 S'o,i[(5y] are the quadratic in SY parts of the action, and the subscripts "0" and "1" denote the action 
respectively without and with one QPS. 

The integrals over fluctuations in eq. (54) can be evaluated exactly only in simple cases. Technically such a 
calculation can be quite complicated even if the saddle point trajectories can be determined explicitly. In our case an 
anal ytic al expression for the QPS trajectory is not even known. Hence, an exact evaluation of the path integrals in 
eq. (p4|) is not possible. 

Below we will present a simple approach which allows to establish the correct expression for the pre-exponent B up 
to an unimportant numerical prefactor. Within our present analysis any attempt to find an explicit value for such a 
prefactor would make little sense simply because the numerical value of A in eq. (^7|) is not known exactly. Also for 
other problems numerical prefactors in the pre-exponent are usually of little interest. Therefore we believe that our 
approach may be useful for various other situations because it allows to establish the correct functional form of the 
pre-exponent practically without any calculation. If needed, with a little extra effort our method may also allow to 
approximately evaluate a numerical coefficient in the pre-exponent. 

In order to calculate the ratio of the path integrals in eq. j54| ) let us introduce the basis in the functional space 
^fc(z) in which the second variation of the action around the instanton S 2 Si[SY] is diagonal. Here the basis functions 
depend on a general vector coordinate z which is simply z — (t,x) in our case. The first N functions are the 
so-called "zero modes" related to the invariance of the instanton action under arbitrary shifts in certain directions in 
the functional space (in our case - shifts of the QPS position along the wire and in imaginary time, i.e. N — 2). Let us 
denote an instanton solution as Y(z). Then the zero mode eigenf unctions are expressed as follows: ^>k(X) = dY/dzk, 
where k < N and the number of zero modes N coincides with the dimension of the vector z. An arbitrary fluctuation 
5Y(z) can be represented in terms of the Fourier expansion 

SY(z) = J2Sz k ^^+ J2 u k * k (z). (55) 

fc=l k k=N+l 

Then we get 

S 2 S [SY] = - J2 A kn u k u* n , 5 2 S 1 [5Y} = - J2 X M 2 , (56) 



k,n=l k=N+l 



where for k < N the Fourier coefficients uu = Sz^ are just the shifts of the instanton position along the k— th axis 
and A/c are the eigenvalues of (5 2 Si [<5y ] . Integrating over the Fourier coefficients one obtains 
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JVSYexp(~6 2 Si[5Y}) 
JVSYexp(-5 2 S Q [SY}) 



dSxi.. / dSxpj 



\ 



, (57) 

(2tt)^ J] A fc 

fc=iV+l 



where is the system size in the k— th dimension. The forimda (|5j) is, of course, not at all new. It just represents 
the standard ratio of determinants with excluded zero modesc3. We will argue now, that with a sufficient accuracy 
in the latter formula one can keep the contribution of only first N eigenvalues. Indeed, the contribution of the 
"fast" eigenmodes (corresponding to frequencies and wave vectors much larger than the inverse instanton size in the 
corresponding dimension) is insensitive to the presence of an instanton. Hence, the corresponding eigenvalues are 
the same for both S 2 Sq and 5 2 Si and just cancel out from eq. ([37|). In addition to the fast modes there are several 
eigenmodes with frequencies (wave vectors) of order of the inverse instanton size. The ratio between the product of 
all such modes for S 2 S\ and the product of eigenvalues for 5 2 Sq with the same numbers is dimensionless and may only 
affect a numerical prefactor which is not interesting for us here. Dropping the contribution of all such eigenvalues one 
gets 



JV5Yexp(-5 2 S 1 {5Y}) 
JV8Yexp{-S 2 S [SY}) 




(58) 



What remains is to estimate the parameters A kk for k < N. For this purpose let us observe that the second variation 
of the action becomes approximately equal to the instanton action, 5 2 Si = \A kk z^ k ~ Sqps, when the shift in the 
A:— th direction becomes equal to the instanton size in the same direction Sz k = ZQ k . Then we find Akk ~ %Sqps / z o k 
and 

det A kn \ k , n<N « J] A kk » (2 ^ QP5)jV . (59) 

n z 2 



Ok 

k=l 

Finally, combining eqs. (|50|), (|54|), ( |58| ) and (|5|) we obtain 

— (n^)(^rf' « 

Here b is an unimportant numerical prefactor. This result demonstrates that the functional dependence of the pre- 
exponent can be determined practically without any calculation. It is sufficient to know just the instanton action, the 
number of the zero modes N and the instanton effective size zofc for_each of these modes. 

Let us also note that a similar observation has already been madeCJ for some local Lagrangians equal to the sum of 
kinetic and potential energies. Here we have shown that the result ( |60| ) holds for arbitrary effective actions, including 
nonlocal ones. Hence, this result can be directly applied to our problem of quantum phase slips in thin superconducting 
wires. In this case we have L\ = 1/T, L2 = X and eq. ( |60| ) yields 

B^^^. (61) 

T0X 

This equation provides an accurate expression for the pre-exponent B up to a numerical factor of order one. As we have 
already discussed, such an accuracy is sufficient for our purposes. We also note that the result (^l|) is parametrically 
different from previous results obtained within a TDGL-type of analysis^ or suggested phenomcnologically in Rcf. [|. 

Finally, it is worth mentioning that we have also compared our eq. (J6^) with the exact results previously obtained 
for various other problems by means of different approaches. Here we will briefly discuss three different examples for 
the sake of illustration. The first example is the problem of a quantum particle in a cosine periodic potential. In this 
case an explicit expression for the inter-well tunneling rate is well known. If we apply our method and evaluate only 
An, we will obtain the tunneling rate which is ~ 40 % smaller than the exact result. If we also evaluate A22 and A2 
and include their ratio into our formula, the result for the tunneling rate will be only 10 % smaller as compared to the 
exact one. This example demonstrates that also a sufficient numerical accuracy in the pre-exponent can be achieved 
without a complicated calculation of the ratio of the determinants. 

Two other examples concern the systems with neslocal in time Lagrangians. Consider, e.g., the problem of quantum 
decay of a particle in the presence of dissipationca. In the limit of strong dissipation this problem was treated by 
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Larkin and OvchinnikovtS who found the exact eigenvalues and, evaluating the ratio of the determinants, obtained 
the prefactor in expression for the decay rate B cx rf^ 2 /m 2 , where 77 is an effective friction constant and m is the 
particle mass. This result would imply that the pre-exponential factor in the decay ratet2l should be very large and 
may even diverge if one,formally sets m — * 0. ■— ■ 

Later it was realizedU that this divergence is artificial. Performing_a. simple one-loop perturbative calculationEil one 
arrives exactly at the same high frequency divergence as in the resulted. It implies that this divergence has nothing to 
do with tunneling, and it is regularized by means of a proper renormalization of the bare parameters in the effective 
action. After that the high frequency contribution to the pre-exponent is eliminated and one findsEJ B oc 1/ ^/rj. [Note 
a misprint in the power of 77 in eq. (8) of Ref. |3l}] This expression does not contain JJjie particle mass m at all. It 
also allowed to fully resolve a discrepancy with the experiments^. Note that the resultEj can also be expressed in the 
form B ~ \fSb/To, where S& oc rj is the instanton (bounce) action and To oc 77 is its typical size. As we have already 
argued, the result in this form can be guessed from eq. (60) without any calculation. [For this particular problem our 
approach allows to even reproduce an exact numerical prefactor.] 

The last example is the problem of Coulomb blockade in normal tunnel junctions in the strong tunneling limit. 
This problem was treated within the instanton technique in Ref. |3^. In this problem each instanton has two zero 
modes which correspond to its shifts in time and fluctuations of its frequency £1. The value of the instanton action is 
well known and the parameters zqi and Z02 for both zero modes can be evaluated directly by means of the approach 
presented above. Each of these parameters is found to depend on one of the zero modes f2. However, the product 
Z01Z02 turns out to be an ^-independent constant. Making use of this fact and integrating over the zero mode 
coordinates in eq. (60), we arrive at the functional form of the pre-exponent derived in eq. (10) of Ref. p$ by means 



of an explicit calculation of the fluctuation determinants, see also Ref. 22 



The above examples demonstrate that our approach allows to easily derive the functional form of the pre-exponent 
in a variety of problems, including those where technically involved calculations appear to be inevitable otherwise. 



V. COMPARISON WITH EXPERIMENT 

Now let us compare our results with experimental findings. Recently Bezryadin, Lau and TinkhamEl reported a 
clear experimental evidence for the existence of quantum phase slips in ultra-thin (with diameters dpsvn to 3 nm) and 
uniform in thickness superconducting wires. Three out of eight samples studied in the experimentsEj showed no sign 
of superconductivity even well below the bulk critical temperature Tq- Furthermore, in the low temperature limit the 
resistance of these samples was found to show a slight upturn with decreasing T. In view of that ope can conjecture 
that these samples may actually become insulating at T — > 0. The resistance of other five samplesEj decreased with 
decreasing T. Also for these five samples no clear superconducting phase transition was observed. 

All three non-superconducting wires (il,i2 and i3) had the normal state resistance below the quantum unit R q , while 
the normal state resistance of the remaining five "superconducting" samples was larger than R q . This observation 
allowed the authorscJ to suggest that a dramatic difference in the behavior of thesej±ttp,group of samples (otherwise 
having similar parameters) caa. be due to the dissipative phase transition (DPT)li3'EIlE3 analogous to that observed 
earlier in Josephson junctional]]. 

Without going into details here, let us just point out that DPT can be observed only provided quantum phase 
slips are easily created inside the wire. The results for 7qps derived in the present paper allow to estimate a typical 
average time within which one QPS event occurs in the sample. Making use of eqs. @, @ and @ we performed 
an estimate of such a time to — 1 /'Jqps for all eight samples studied in Ref. ^0|. In this experiment the samples were 
fabricated from MoygGe2i alloy. For our estimates we will use the value of the density of states No = 1.86 x 10 13 
sec/m 3 for clean Mo, which can be extracted from the specific heat data. The resistivity of the material was measured 
to be p — 1.8 /.tfi/m, the superconducting critical temperatese is Tc — 5.5 K. With these numbers we obtain the 
coherence length £ ~ 7 nm in agreement with the estimatec! The results for to are summarized in the following 
Table: 



sample 


R/d, kO/nm 


So 


£oU=i, sec 


ioU=2, sec 


il 


0.122 


7.8 


HT 11 


10- 8 


i2 


0.110 


8.7 


io- n 


10~ H 


i3 


0.079 


12.7 


10" y 


10- 4 


si 


0.038 


25.1 


10- 4 


10 B 


s2 


0.028 


33.7 


1 


10 i4 


s3 


0.039 


22.6 




10 b 


ssl 


0.054 


15.4 




HT 1 


ss2 


0.044 


19.6 


io- B 


10 2 
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The action S is denned by Eq. @ with A=l. The typical QPS time t 



is very sensitive to the particular value of the factor A, therefore here we present two estimates corresponding to 
A = 1 and A = 2. 

In spite of remaining uncertainty in the prefactors some important conclusions can be drawn already from the 
above estimates. For instance, we observe that for both A = 1 and A — 2 the QPS rate 7qps = 1/to is very high 
(as compared, e.g., to the typical experimental time scale ~ 1 sec) in the "insulating" wires il, i2 and i3. This fact is 
fully consistent with the observationsEHl: numerous quantum phase slips occurring in these wires completely destroy 
the phase coherence and, hence, superconductivity is washed out. Thus, non-superconducting behavior of these three 
samples should be due to quantum phase slips. ■— ■ 

On the other hand, the QPS rate is notably lower for all the "superconducting" wired 2 !! Possible interpretation 
of the experimental results for the samples sl-ss2 depends strongly on the value of A. E.g. for A = 1 the QPS rate 
is high enough practically-|in all samples. In this case quantum phase slips should in principle-be important also for 
"superconducting" wirescl Then one can indeed relate the behavior of these samples to DPTlLa, as a result of which 
quantum phase slips are bound in pairs and, hence, quantum fluctuations are strongly suppressed. 

If, however, one chooses A = 2 the QPS time for the samples sl-ss2 turns out to be very long, much longer than 
the experimental time. Then the QPS effects should be irrelevant, and one would expect these samples to show 
a superconducting behavior, perhaps with the_renormalized critical temperatures. This conclusion would also be 
consistent with the experimental observations 2 !]. 

Finally, let us note that all the above estimates are performed in the limit T = 0. This is correct if temperature is 
considerably below Tq. Otherwise the expression for the QPS action Sqps needs to be modified. 

In conclusion, we have developed a detailed microscopic theory of quantum phase slips in ultra-thin homogeneous 
superconducting wires. We hapj derived the effective QPS rate for such wires and evaluated this rate,for the systems 
studied in recent experimentsEj. Our results are fully consistent with the experimental findingscJ which provide 
perhaps the first unambiguous evidence for QPS in mesoscopic metallic wires. 



APPENDIX A: 



Let us collect some rigorous expressions for the "susceptibilities" xe, X.J-, Xl and xa- In Ref. [l6| these quantities 
have been related to the so-called polarization bubbles fo^-fo and ho- In the interesting for us diffusive limit Al/vp <C 1 
these polarization bubbles are defined by the equationsEEl 

f rf 3 k A 2 

/„=rV - — ^(F(w + w v ,a + k)F(u} v .k)),, = ttAoTV ; (Al) 

^ J (2tt) 3 x v ' H ; K ' " dls WW (W + W + Dq 2 ) K ' 

go=T y [ ^(G( W + ^ )q + k)GK.k)), =~N + nN Q TY WW' - u(w + u v ) 

^ J (2tt) 3 x v ; v //dls ^ WW'(W + W' + Dq 2 ) v ' 



h tV f d3k Ini _l _i_ i M\ wtV WW' + + u v ) 

h a = T^ J -j—^ {G(tJ + uj v ,q + k)G(uj v ,k)) d . s = -ttNqT^ 



WW(W + W + Dq 2 
In these equations we use the notations 



(A3) 



W = y/u}* + A 2 , W' = \J (uj y + u) 2 + A 2 , (A4) 

w v = -kT(2v +1) and u = 2itTn, where v and n are arbitrary integers, and (. . .)di s implies disorder averaging. 
Let us start from the function xe defined as 

2e 2 

Xe = a-(/o+5o). (A5) 

q 2 

At T = we find 
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XE = < 



+00 



A 2V1+Z 2 jj" ^■ v /(t 2 -l) 2 +4t 2 /(l+x 5 )(l+t 2 +^/(t 2 -l) 2 +4t 2 /(l+x 2 )) 



l-t 2 + ^/(t 2 -l) 2 +4t 2 /(l+:r 2 ) 



y = o, 



A 2j 4y 



, a; = 0, y > 1, 



i(^-^ + ^farctan( v &)), * = 0, y < 1. 



Here and below we set x = uj/2A and y = Dq 2 /2A. In the limit |w| » 2A we obtain 

<7 



X£ - 



|w| +Dq 2 



At small w and q we get 



Xb = 



8A 



3 / w \ 2 8 £>g 
~ 8 V2AV ~3^~2A 



2 1 



It is also possible to evaluate X£ at uj — 0, q = and arbitrary T: 

A A 



Xs(0,0) 



tanh 



8A ^ 2T 2T(cosh 



2TV 



For A < T we find 



Xb(0,0) = 



ttctA 2 



96T 3 ■ 

In this limit also a more general expression for arbitrary lu and g can be obtained: 



Xe 



4aA 2 Dq 2 



\uj\+Dq 2 (L0 2 -D 2 q 4 ) 2 
2aA 2 (co 2 + D 2 q 4 ) 



1 , \uj\ + Dq 2 



M(w 2 -£>V) 2 



2 2ttT 



4ttT 
1 
2 



1 



0A 2 , , 

2ttT(u; 2 - £>V) ^2 + 



27rry 



where ^(x) is the digamma function. 

The remaining x-functions are evaluated analogously. Consider the function 



Xj = 4e 2 / . 



After straightforward algebra we obtain 



XJ = < 



2 ln(|x|+vW ) 

4e 2 iv n f — - l <y + V^) 



4^0 £ 



: arctan 



' ~ 3 \2A/ ~ 4~2A~ 



In the limit of low frequencies and wave vectors one has 

XJ = 2e 2 iV 
At uj = 0, q = and for arbitrary T we get 

X J ^2 7 r e 2 iV A 2 Tg ( ^ + 1 A2)3/2 ={ 

where ((3) ~ 1.202. In the limit T > A and u> ^ wc find 



= 0, 
x = 0, y > 1, 



a; = 0, y < 1. 



2e 2 N , T < A, 

7C(3) e 2 jV A 2 T a 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 
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